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We discuss the topological p x + ip y superfluid phase in a 2D gas of single-component fermionic 
polar molecules dressed by a circularly polarized microwave field. This phase emerges because the 
molecules may interact with each other via a potential Vo(r) that has an attractive dipole-dipole 
1/r 3 tail, which provides p-wave superfluid pairing at fairly high temperatures. We calculate the 
amplitude of elastic p-wave scattering in the potential Vo(r) taking into account both the anomalous 
scattering due to the dipole-dipole tail and the short-range contribution. This amplitude is then 
used for the analytical and numerical solution of the renormalized BCS gap equation which includes 
the second order Gor'kov-Melik-Barkhudarov corrections and the correction related to the effective 
mass of the quasiparticles. We find that the critical temperature T c can be varied within a few 
orders of magnitude by modifying the short-range part of the potential Vo(r). The decay of the 
system via collisional relaxation of molecules to dressed states with lower energies is rather slow due 
to the necessity of a large momentum transfer. The presence of a constant transverse electric field 
reduces the inelastic rate, and the lifetime of the system can be of the order of seconds even at 2D 
densities ~ 10 9 cm -2 . This leads to T c of up to a few tens of nanokelvins and makes it realistic to 
obtain the topological p x + ip y phase in experiments with ultracold polar molecules. 



I. INTRODUCTION 



The recent breakthrough in creating ultracold diatomic 
polar molecules in the ground ro-vibrational state [TJ [5] 
and cooling them towards quantum degeneracy pQ has 
opened fascinating prospects for the observation of novel 
quantum phases [3HS]- A serious problem in this direc- 
tion is related to ultracold chemical reactions, such as 
KRb+KRb=>K2+Rb2 observed in the JILA experiments 
with KRb molecules [10] , which places severe limitations 
on the achievable density in three-dimensional samples. 
In order to suppress chemical reactions and perform evap- 
orative cooling it has been proposed to induce a strong 
dipole-dipole repulsion between the molecules by confin- 
ing them to a (quasi)two-dimcnsional (2D) geometry and 
orienting their dipole moments (by a strong electric field) 
perpendicularly to the plane of the 2D translational mo- 
tion [TT]. Nevertheless, in order to prevent the approach 
of colliding molecules at short separations, where the 
chemical reactions occur, and to proceed with evapora- 
tive cooling, one has to have a fairly strong confinement 
to the 2D regime [T2TH3] . The suppression of chemical re- 
actions by nearly two orders of magnitude in the quasi2D 
geometry has been demonstrated in the recent JILA ex- 
periment (15] . It should be noted, however, that these re- 
actions are not expected to occur for all polar molecules 
of alkali atoms [IB], on which experimental efforts are 
presently focused. While present for KRb molecules and 
for the molecules containing a Li atom, they are ener- 
getically unfavorable, for example, in the case of NaK 
molecules. We thus expect that future investigations of 
many-body physics will deal cither with molecules which 
do not undergo ultracold chemical reactions, or otherwise 
are strongly confined to the 2D regime. 



One of the challenging goals in the studies of many- 
body physics is the creation of the topological p x + ip y 
superfluid phase for identical fcrmions in two dimensions 
(2D) (see [T7] for a review). This phase, first discussed 
in relation to superfluid 3 He and the fractional quantum 
Hall effect [TS1 HH] , has exotic topological properties at 
positive chemical potential [i > 0, i.e. in the BCS regime. 
The topological nature of this superfluid phase supports 
a gapless Majorana mode at the boundary to a vacuum, 
and the quantized vortices in the superfluid carry local 
zero energy Majorana modes on their cores [Ti?l [213] , These 
Majorana modes are predicted to cause the vortices to 
obey non-abelian exchange statistics, which has poten- 
tial applications in topologically protected quantum in- 
formation processing |21) . There is significant interest 
in finding physical realizations of this topological super- 
fluid phase in which the exotic physics of these Majorana 
modes can be detected. 

The p x + ip y topological phase can be the ground state 
of ultracold identical fermionic atoms interacting via a 
short-range potential [17] , However, away from a p- wave 
Feshbach resonance the superfluid transition tempera- 
ture T c is vanishingly low. On approach to the reso- 
nance it increases but the system becomes unstable due 
to the formation of long-lived diatomic quasibound states 
and their collisional relaxation into deep molecular states 
[2"2l |2"5] . A stable p x + ip y state has been recently pro- 
posed for fermionic polar molecules with a large dipole 
moment [5]. When these molecules are confined to a 2D 
geometry and dressed by a microwave field (MW) which 
is nearly resonant with the transition between the low- 
est and the first excited rotational levels, they acquire an 
attractive 1/r 3 dipole-dipole interaction. This leads to 
superfluid pairing of p x + ip y symmetry, and due to the 
anomalous contribution of the 1 /r 3 tail to the scattering 
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amplitude, already in the BCS limit the superfluid tran- 
sition temperature can be made a sizeable fraction of the 
Fermi energy. At the same time, collisional decay pro- 
cesses remain sufficiently slow to allow the experimental 
realization of this phase |SJ . Several other ways in which 
a stable p-wave coupled superfluid phase may be obtained 
have been proposed in recent years |24H30j . 

In this paper we present a complete analysis of the BCS 
limit of the p x + ip y phase of MW-dressed fermionic polar 
molecules. In particular, using the many-body perturba- 
tive approach up to second order, we derive a relation 
for the transition temperature T c . It is then shown how 
T c and the collisional stability may be manipulated by 
tuning the short-range part of a MW-induced effective 
molecule-molecule interaction potential or by applying 
a constant electric field perpendicular to the plane of 
the translational motion. Our analysis is confirmed by 
numerical calculations, which are also extended to the 
regime of moderately strong interactions. 

The paper is organized as follows. In Section II we 
discuss a single polar molecule in the presence of a cir- 
cularly polarized nearly resonant microwave field and a 
constant electric field as shown in Fig. 1. We then show 
that two microwave-dressed polar molecules undergoing 
a 2D translational motion may interact with each other 
via a potential Vo(r) which has a repulsive core, a po- 
tential well, and an attractive 1/r 3 tail. In Section III 
we find the amplitude of elastic p-wave scattering of par- 
ticles in the potential Vo(r), which takes into account 
both the anomalous scattering due to the 1/r 3 tail and 
the short-range contribution. Section IV is dedicated to 
the analysis of the decay of the gas due to collisional re- 
laxation of the molecules to lower dressed states. The 
presence of a constant electric field is found to reduce 
the inelastic rate, and the lifetime of the system can be 
of the order of seconds even at 2D densities ~ 10 9 cm~ 2 . 
In Section V we present the renormalized BCS gap equa- 
tion which takes into account the second order Gor'kov- 
Melik-Barkhudarov processes and the effective mass of 
the quasiparticles. In Section VI we obtain analytical 
and numerical solutions of the gap equation, show that 
the ground state has p x + ip y symmetry, and reveal the 
influence of the short-range part of the scattering poten- 
tial Vo(r) on the critical temperature. In Section VII we 
conclude, emphasizing that it is realistic to obtain the 
topological p x + ip y phase for microwave-dressed polar 
molecules in the 2D geometry, with a critical tempera- 
ture of up to a few tens of nanokelvin. 



II. INTERACTION POTENTIAL FOR 
MW-DRESSED POLAR MOLECULES 

The microwave dressing of polar molecules has been 
proposed to tune the molecule-molecule interaction po- 
tential [31 and to form a repulsive shield for suppressing 
inelastic losses |32]. The rotational part of the Hamilto- 
nian of a single polar molecule in a circularly polarized 




FIG. 1: (color online). Polar molecule (brown and green cir- 
cles) in the considered configuration of fields. The ac MW 
field rotates with frequency cj in the plane orthogonal to the 
dc field 



MW-field E ac and a constant electric field Edc reads: 

i/ = BJ 2 -d-(E dc + E ac ), (1) 

where J is the operator of the rotational moment, B is the 
rotational constant of the molecule, and d is the operator 
of the dipole moment. In the following we assume that 
Edc is parallel to the z-axis, and the vector E ac rotates 
in the {x,y} plane (see Fig. 1). 

In the absence of electric fields, good quantum num- 
bers are the rotational moment J and its projection 
on the z-axis, Mj, and the diagonal matrix element 
(J, Mj\d\J, Mj) is equal to zero. The cfc-field preserves 
Mj but couples states with different J. We will consider 
a rather weak cfc-field so that the parameter 



dE dc 
B 



(2) 



is small, where d is the molecular dipole moment. Within 
the manifold that includes the original states with J = 
and J = 1, the new basic states will be called | ^oo) ? 

|$io), and |$n). The transition matrix elements 
of the operator of the dipole moment, which are needed 
for further discussions, are given by [31] : 



d c = |($oo|<M$i±i 
d g = K$oo|4|$oo)| 



d/3 

T' 



de = K$l±l|d,|*l ±1 ) 



d/3 
10' 



(4) 
(5) 



where we have omitted terms that are higher order in /3. 

Consider the application of a circularly polarized MW- 
field with frequency u) close to the frequency luq of the 
transition between the states |$oo) and | <&xi) - For the 
Rabi frequency £Ir = d c E ac /h and the frequency detun- 
ing 5 = uj — luo satisfying the inequality \S\, flu <C wq, the 
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rotating wave approximation is valid and the MW-field 
couples only the states |$oo) an d |$n)- The resulting 
eigenstates are: 



interaction between the molecules: 



a|$ 00 ) +6|$u), 
&|$oo) - a|$n), 



(6) 
(7) 



with a 



~W R ,b 



A 2 , and A 



(5+^/S 2 + 4fl R ) /2. We will consider (5 > O fl , and choose 
<5 > such that the energy of the state | +) lies above the 
energies of the states |— ) and |$i_i). Ramping the MW 
field on adiabatically then ensures that the ground state 
|0, 0) evolves into the state |+). Thus all molecules can 
be prepared in this state. 

We remark that it was demonstrated in Ref. [33] that 
polar molecules may be prepared in any hyperfine sub- 
state of the ground ro-vibrational state. In particular, we 
shall assume that the molecules are prepared in the state 
with maximum projection of the magnetic moment such 
that the circularly polarized MW field does not cause 
mixing of hyperfine levels. 

We now consider two MW-dressed polar molecules un- 
dergoing a 2D motion in the {x, y} plane and separated 
by a distance r which greatly exceeds the radius of the 
ordinary van der Waals interaction potential. The pres- 
ence of the dc and MW-fields introduces a dipole-dipole 



H, = 



di-d 2 -3(di-f)(d 2 -f) 



(8) 



Assuming that the molecules are at a fixed separation 
r we employ the Born-Oppenheimer approximation and 
determine an effective intermolecular potential induced 
by the dc and MW-fields. For this purpose we have 
to diagonalize the Hamiltonian H = H + Hd- The 
dipole-dipole interaction couples the molecular states 
|s) = |$ 00 ), |e) = |$u), and |e) = |$i_i), whereas the 
state |$io) remains decoupled from them. The wavefunc- 
tion of a two-molecule dressed state can be cither sym- 
metric (even parity) or antisymmetric (odd parity) with 
respect to permutation of the molecules. For studying 
topological superfluids we will be interested in the states 
of even parity. We then have a basis of five two-particle 
states: 

(|e, g ) + | g ,e)) (\e,g) + \g,e)) (\e, e) + \e, e)) 

\9,9), 7^ j 7E j |e,e), - 



V2 



V2 



V2 



(the state |e, e) is decoupled from these states). The 
Hamiltonian acting on the spinor of these five states is 
represented by a matrix 32 : 



H = ns 



( (3 2 /3x 3 y/2Q R /S 

V2n R /S (/3 2 /10 - 1/2) /x 3 - 1 3/2x 3 

3/2x 3 (/? 2 /10- 1/2) /x 3 

V2Q R /5 

Q, R /5 





V2~n R /5 

1 

3/3 2 /100x 3 - 2 




\ 


n R /5 
o 

3/3 2 /100x 3 - 2 / 



where x = r/r§, r$ = (d^/HS) 1 ^ 3 , and we omit terms 
that are higher order in ft. Diagonalizing the matrix 
([9]) we obtain a number of effective intermolecular po- 
tentials (potential curves) V e!! {r) for MW-dressed polar 
molecules, depending on the states of the molecules at 
an infinite separation (see Fig. 2). 

For two molecules which at r — ¥ 00 are in the state 
|+), the long-range tail of such a potential denoted here 
as Vo(r) is 



V (r -> 00) = - 

with the length scale r* given by 
md 2 (n R /6) 2 



h 2 



m r 



3 ' 



(10) 



m 2 i + A{n R /8f 



{l -3^ 2 x 



49/2160 



'7+13 v /l + 4(^/(5) 2 

6on R /s 



}< 
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and decreasing with an increase in /3. For small £Ir/S and 
sufficiently large j3 the parameter becomes negative, 
which leads to a repulsive tail of Vo(r). In this case the 
de-field dominates over the effect of the MW-dressing, so 
that on average the molecules are perpendicular to the 
{x, y} plane and exhibit a repulsive dipole-dipole interac- 
tion at large r. In the following we consider only the case 
where is positive and the long-range tail of Vq (r) is at- 
tractive. Then the quantity is a measure of the radius 
of the centrifugal barrier experienced by the (fermionic) 
molecules. 

At smaller separations, the dipolar interactions be- 
tween the molecules cause them to depart from the state 
|+). This occurs when the characteristic interaction en- 
ergy d 2 /r 3 becomes larger than the detuning h\S\, setting 
the lengthscale r$ = (dl/hS) 1 ^ 3 entering the matrix 
We assume that r$ is larger than the confinement length 
in the z-direction, l z , so that the interaction is 2D. The 
potential curves of even parity are shown in Fig. 2. The 
potential Vo(r), being attractive at r > r& with a long- 



4 



0.01 



0.005 - 



us; 




-0.005 



FIG. 2: (color online). Potential energy curve Vo(r) for two 
|+) state molecules, computed for fijj = 0.255. The black 
(solid) curve corresponds to the value of the perpendicular 
dc field P = 0, the blue (dashed) curve to /? = 0.1, and the 
red (dash-dotted) curve to f3 = 0.2. Anticrossings of Vo(r) 
with other field-dressed levels of even parity occur at distances 
r ~ r$, as shown in the inset for Qr = 0.255 and /3 = 0. 



range 1/r 3 tail, has a repulsive core for 



< 



The 



repulsive core prevents low-energy molecules from ap- 
proaching each other at distances r < r$ and suppresses 
inelastic collisions, including ultracold chemical reactions 
observed at JILA in experiments with KRb molecules 
[ID] . Note that in this respect the MW-dressing of po- 
lar molecules in the 2D geometry can be used for their 
evaporative cooling. 

Actually, the potential curves of even parity, in partic- 
ular Vo(r), may cross with potentials of odd parity. How- 
ever, inelastic transitions between the states of different 
parity are not possible in two-body collisions, and three- 
body collisions accompanied by such transitions will be 
suppressed at the low densities that we consider. 



III. ELASTIC SCATTERING 

We now consider elastic scattering of fcrmionic 
molecules (each in the state |+)) undergoing 2D transla- 
tional motion and interacting with each other via the po- 
tential Vo(r), with the attractive dipole-dipole tail (10 1. 



At ultralow energies the leading scattering is in the p- 
wave scattering channel. For the investigation of super- 
fluid pairing we need to know the off-shell scattering am- 
plitude defined as 

/(k', k) = / expHk' • r)V (r)Mr)d 2 r, (12) 



where ipk(r) is the true wavefunction of the relative mo- 
tion with momentum k. The p-wave part of /(k', k) is 
f(k', k) expi0, where <j) is the angle between the vectors 



k and k', and the partial p-wave amplitude is given by 

poo 

f(k',k)= J 1 (k , r)V (r)4>(k,r)2Trrdr, (13) 
Jo 

with Jx being the Bessel function. The wavefunction of 
the relative p-wave motion, 0(fc,r), is governed by the 
Schrodinger equation 



d 2 ip 1 dip -0 



ni 



dr 2 



r dr 



h 2 k 2 T 



V>. (14) 



The full on-shell scattering amplitude is obtained from 
equation (12) at |k'| = |k|, and its p-wave part follows 
from Eq. (13) at k' = k, f(k) = f(k, k). It is related to 



the p-wave scattering phase shift 8[k) by [34] 

„.,. 2h 2 , , , 4h 2 tanS(k) , . 

/(*) = r exp(2W (k )-l = { ' (15) 

mi m 1 — itano(fc) 



In Eq. ( 13 ) the wavefunction tp(k, r) is normalized such 
that for r — > oo we have 

^(k,r) = J 1 {kr)- l ^-H 1 {kr), 



where Hx = J\ + iN\ is the Hankel function, and N\ is 
the Neumann function. It is, however, more convenient 
to normalize the radial wavefunction in such a way that 
it is real and for r — > oo one has 

V>(fc,r) = [ Ji (fcr)-tan 5{k)Nt (kr)} oc cos(/cr-37r/4+(5) . (16) 

One easily checks that tp(k, r) = tp(k, r)/(l — i tan S(k)). 
Using this relation in Eq. (13) we see that the off-shell 



scattering amplitude can be written in the form: 

f(k', k) 



f(k',k) = 



1 — i tan S(k) ' 



(17) 



where f{k', k) is real and is given by Eq. ( 13 ) with ip(k, r) 
replaced by tp(k,r). For k' — k one has 



f(k,k) = f(k) = -(4h 2 /m)tanS(k). 



(18) 



To find the p-wave scattering amplitude in the limit 
fcr* C 1 we divide the range of distances into two parts: 
r < r and r > r , where r lies in the interval r + <C 
ro <C fc _1 (see Fig. 3). As we will see, this separation 
reflects the presence of two contributions to the scattering 
amplitude, the short-range contribution (r % r*) and the 
anomalous contribution (r ~ k^ 1 ). 

In region I where r < ro, the p-wave relative motion of 
two fermionic molecules is governed by the Schrodinger 
equation with zero kinetic energy: 



h 2 (d 2 i\)i l#j V/ 



dr 2 



+ 



r dr 



+ V ipi = 0. 



(19) 
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- 










X ~-l/r 3 . 



For the asymptotic representation of the relative wave- 
function, chosen in Eq. (16 1, we have 



G(r,r') 



m f ^(r^N^kr), r > r' 
4h 2 \ t/>ff(r)Ni(kr'), r<r' 



(25) 



Then, substituting the Green function (25 1 into equation 



(231 and taking the limit r — > oo we obtain: 



tan<5 (1) (fc)=tan < 5 / (fc)- 



4h' 



<Pn( r ) V (r)2Trrdr.{26) 



FIG. 3: (color online). The division of ranges into regions / 
(r < ro) and II (r > ro). The length scales r t and ro 2> r* 
are indicated by a circle and a dashed line respectively. 



Using Eq. (22 I for tan 5/ and calculating the integral in 



Eq. ( 26 ) we find: 



w 3 8 



In— +2C--+nA 
r 2 



> (27) 



At distances where the interaction potentia l is behaving 
as Vo(r) — ~h 2 r^,/mr 3 , the solution of Eq. (19) reads: 



■0/ (r) oc 



AJ 2 2J- 



N 9 . [ 2J — 
r 



(20) 



where the constant A is determined by the behavior of 
Vo(r) at shorter distances. 

In region II, at r > ro, the relative motion is prac- 
tically free and the potential Vo(r) can be considered a 
perturbation. To zero order we then have for the relative 
wavefunction: 



^(r) = Ji(fcr) -tsm5i(k)Ni(kr), 



(21) 



and we omitted terms that contain higher powers of k. 

Adding the second order contribution we have for the 
relative wavefunction: 



ipy/(r)= / G(r,r')V {r')2Trr'dr' 



x / G(r',r")V (r")i}Y/{r")2Trr"dr". (28) 

J To 



Then, using Eq. ( 25 1 and taking the limit r — > oo we 
see that including the second order contribution, the scat- 
tering phase shift becomes 



where the scattering phase shift Si is due to the interac- 
tion between particles in region I. Matching the logarith- 
mic derivatives of ipi and ifrffl at r — ro we obtain 



tan Si = - 



l-^ hC-^ + nA- In * 
ro V 2 r* 



(22) 



where C = 0.5772 is the Euler constant, and we took into 
account that ro 3> r* and kro <C 1. 

We now include perturbatively the contribution to the 
p-wave scattering phase shift from distances r > ro (re- 
gion II). In this region, to first order in Vo(r) the relative 
wavefunction is given by 

/>oo 

^n\r)=^n\r)- G(r, r')Vo{r')^ {r')2Ttr' dr\ (23) 



where the Green function for the free p-wave motion 
obeys the radial equation: 



tan£(fc) = tan5 (1) (/c) j 



tPnir) V (r)2TTrdr 



x / N 1 {kr')Vo{r')iljf I ) (r')2TTr'dr' . (29) 



As we are not interested in terms that are proportional 
to k 3 or higher powers of k and in terms that behave as 
(&r*) 2 r*/ro, we may omit the term tan Si( k)N i(kr) in 



the expression for ipjj(r). Then equation (291 reduces 
to 



tan^(fc) = tan<5«(fc) - f°° iM dx x 



9 / t2 

^x{N (x) J 2 (x) - N^J^x)) - ^N {x)J!(x) 



+ In 1 (x)J 2 (x) - — 

O TTX 



(30) 



m \ ar z r ar r z I lur 



For the first four terms in square brackets we may put 
the lower limit of integration equal to zero. Then, using 
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the relations: 



IV. INELASTIC COLLISIONAL PROCESSES 



'j*( x )N 1 {x)^ = -± 

X 4-7T 

dr 1 

J?ix)J 2 (x)N (x)— = ± 

q X 07T 



jf(x)N (x) 



dx 1 
a; 2 16ir' 

J 2 {x)J 2 {x)N l {x)% = - 



16tt' 



frrn 



t2 , N cte 1 C In 2 1 , , 



we obtain 

tan (?(*)= tan^ 1 ^-^? 



7 
12 



-C-ln2+ln fcr 



Substituting taneJ^^fc) from Eq. (27) we eventually ar- 
rive at the scattering phase shift 

tan 5 (k) = -kr* — lnpkr*, (31) 

o 8 

where p = exp{3C - In 2 - 11/12 + tvA} ~ 1.13exp(7rA). 
Using equations ( 15 1 and ( 18 ) we then immediately ob- 



tain the on-shell scattering amplitude f(k) and the am- 
plitude f(k). Note that f(k) is conveniently represented 
as a sum of two terms: f(k) — fi(k) + /2(h), where 



fi(k) 



3 m 
--.2 



7T ft 

2 m 



(32) 
(33) 



so that in the low-momentum limit the term fi (k) is dom- 
inant. The related tangent of the scattering phase shift 
is contained in the first order contribution from distances 
r > ro , given by the second term on the right-hand side of 
Eq. (26) in which one keeps only Ji(fcr) in the expression 



for ipii(r). This says that f\{k) — — (4ft /m)tsxi5{k) 



J r J 2 (kr)Vo(r)2Trrdr and for Vq given by equation (10) 



the amplitude fi(k) is reduced exactly to the result of 
Eq. ph. 



The off-shell scattering amplitude f(k',k) defined by 
Eq. (l7] ) can also be written as fi(k', k) + f2(k' , k), and 
the leading low-momentum contribution is given by 



fi(k',k) 



Ji(A'r)J l (Jfer)V r o(r)27rrdr 



irh z 



kr*F 



1 1 fc 2 

2' 2' ' fc' 2 



(34) 



where we took into account that for r > r$ the interac- 
tion potential has the form Vo(r) = —h 2 r*/mr 3 . The 
quantity F in Eq. (34) is the hypergeometric function, 
and the result is written for k < k'. For k > k' one 
should interchange k and kl . 



For the considered case of 5 > the dressed molec- 
ular state |+) is higher in energy than the states |— ) 

and I <J> 1 1 ) . Therefore, molecules in the dressed state 

|+) may undergo pair inelastic collisions in which one 

or both are transferred to the state |— ) or | <3> x 1 ) - For 

^ <5) the released kinetic energy is ~ TiS and it can 
cause both molecules to escape from the sample. The 
kinetic energy release requires a momentum transfer of 
~ h/Xs with \$ = yjfb/m5. For \$/rs ^ 1 the particles 
cannot approach each other sufficiently closely to allow 
the required momentum exchange, and one anticipates a 
reduction in the loss rate. The same condition can be de- 
rived semiclassically as the condition of adiabatic motion 
in the potential Vo(r). 

In order to go beyond this limit and deter- 
mine the loss rate for A 5 approaching rg, we have 
solved the full two-body scattering problem, involv- 
ing states of even parity which at an infinite separa- 
tion are: (|+), |+», (|+), |-», (|+>, |$i_i)), (|->, |-», 
and (|-),|$i_i)) [the state (|$i_i), is decou- 

pled]. We calculate numerically the probabilities Pi 
that two |+)-state molecules with relative angular mo- 
mentum I are scattered into any outgoing channel in 
which at least one of them is in the state |— ) or 
1 1,-1). This corresponds to non-adiabatic transitions 
from the potential Vo(r) to the other potentials shown 
in the inset of Fig. 2. The Hamiltonian term which 
causes these non-adiabatic transitions is the kinetic en- 
ergy (Laplacian) term. Defining the spinor Xi( r ) — 
y/ripi(r), where the index i labels the two-particles 
states (|+),|+», (|+>,|-», (|+>,|*i-i», (|-),|-}),and 
(|— ), |$i_i)), we obtain that the Laplacian term acts on 
Xi(r) as 

Lx(r) = Lij{r)xj{r). 
The matrix Lij{r) is diagonal, Lij(r) — Li5ij : with 



Li(r) 





f d 2 


I 2 


-l/4\ 


m 


\ dr 2 


h — 


r 2 J 


* 


f d 2 




-2) 2 - 


m 


\ dr 2 







i = 1,2,4, (35) 



i = 3,5.(36) 



We note that coupling between these scattering channels 
will affect also the elastic scattering of two molecules in- 
cident in the i — 1 channel, in such a way that the scat- 
tering amplitude will differ for molecules in the I = +1 
and I = — 1 channels. This difference is a consequence 
of the fact that time reversal symmetry is broken by the 
circularly polarized MW field. The change in the scat- 
tering amplitude is small for rg <^ r* , which is valid in 
the situations studied below, so it will be a small effect. 
Nevertheless, we note that it will cause a slight differ- 
ence in the energetics (and the transition temperatures) 
of the p x + ip y and p x — ip y phases, providing a small 
symmetry breaking perturbation that will favour one of 
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FIG. 4: Inelastic rate constant a as a function of kr, for 
Q.r = 0.255, rs/Xs = 8, and /3 = 0. 

these two phases. Which one is favoured depends on the 
handedness of the circularly polarized MW field. 

Taking into account that two molecules are lost in each 
inelastic collision, and writing the molecule loss rate as 

n = —an , (37) 

for the 2D inelastic rate constant we obtain: 

a = 4h/mJ2 p i- ( 38 ) 
I 

This rate constant at zero static electric field Ed c has 
been calculated in Ref. 8J. As in 8], we treat parti- 
cles of the outgoing channel as 2D. This is surely valid 
for \$ ^> l z , where l z is the size of the molecule wave- 
function in the tightly confined z-direction (amplitude of 
zero point oscillations). In this case the energy release 
is insufficient to allow transitions to excited states in the 
z-dircction. On the other hand, the 2D treatment of out- 
going particles is relevant in the opposite limiting case of 
Xs -C l z where they have a very small angle ~ /l z out 
of the 2D plane. Thus, this approach should give a good 
result also for the intermediate case Xs ~ l z . 

The dependence of a on rs/Xs at a constant kr,, ob- 
tained in Ref. [8] for (3 — 0, shows the general trend 
of a reduction of inelastic losses with increasing rs/Xs, 
which is consistent with the semiclassical expectations. 
In addition, there is a dramatic modulation of the in- 
elastic scattering rate, arising from an interference of in- 
coming and outgoing waves in the scattering potential. 
We recover these results and investigate the dependence 
of a on kr,. Fig. 4 shows a as a function of kr, for 
Qr/S — 0.25, rs/Xs = 8, and j3 = 0. The dependence 
a oc (kr,) 2 expected in the low-momentum limit is ob- 
served for fcr* < 0.1. 

We also consider the influence of a static electric field 
on the inelastic losses. In Fig. 5 we present the inelastic 
rate constant a versus the ratio rs/Xs for various values 
of (3 at kr, = 1. One sees that an increase in j3 shifts 




i_. i , i i i , ... lL_. i , i , , J . ^ 

5 10 15 20 25 



FIG. 5: (color online). Inelastic rate constant a as a func- 
tion of rs/\s for flft — 0.255, kr, — 1, and = (black, 
solid), /3 = 0.1 (blue, dashed), and j3 — 0.2 (red, dash-dotted). 
The critical value of rs / Xs above which a bound state of two 
molecules exists is indicated with corresponding vertical lines. 



the interference minimum of a towards larger values of 
rs/Xs and leads to lower values of a at these points. For 
example, for /3 — 0.2 the (first) interference minimum is 
located at rs/Xs = 15.3 with a ~ 4 x 10~ 5 h/m, whereas 
for /3 — the minimum is at rs/Xs — 10.7 and the corre- 
sponding inelastic rate constant is a ~ 3 x 10~ 4 ~h/m. 

The characteristic lifetime of the sample is r ~ (an)" 1 , 
and for a deeply degenerate molecular Fermi gas one 
should take the rate constant a at k equal to the Fermi 
momentum kp- We then see that even for fc^r, = 1 the 
lifetime r of, for example, Na 40 K molecules ranges from 
20 to 2 s when increasing the density from 10 8 to 10 9 
cm~ 2 and considering the first interference minimum at 
= 0.2 (see Fig. 5). 

We should avoid the presence of bound states of two 
molecules in the potential Vo(r). Otherwise three-body 
recombination will lead to a rapid decay of the gas [8]. 
A dimensional estimate for the three-body recombina- 
tion rate constant gives a Iec ~ (Tir* 2 / 'm)(fci?r*) 4 . The 
corresponding decay time is r rcc ~ (a rec n ) , and it 
can be as short as milliseconds or even lower at densities 
n ~ 10 8 — 10 9 cm~ 2 for not very small values of kpr*. As 
we will see below, an efficient superfluid pairing and suf- 
ficiently high superfluid transition temperature require 
kpr, approaching unity. 

Bound molecule- molecule states in the potential Vb(r) 
appear under an increase in rs/Xs- Critical values of 
rs/Xs below which bound states are absent for a given 
(3, are shown by dashed vertical lines in Fig. 5. Thus, 
to the left of these lines the three-body recombination is 
absent. This is in particular the case at the first inter- 
ference minimum for which the lifetime of the gas due to 
two-body inelastic collisions was estimated above. 



V. SUPERFLUID PAIRING AND THE GAP 
EQUATION 

We now discuss superfiuid pairing in the 2D gas of 
MW-dressed fermionic polar molecules in the internal 
state |+), interacting with each other via the potential 
Vo(r). The Hamiltonian of the system reads: 



d 2 rft(r) 



2m 



\ [ d 2 rd 2 r'^(r)^(r')V a (r^r')^(r')'tP{r), (39) 



u 



where $(r) is the field operator of the dressed molecules, 
and n is the chemical potential. The pairing is due to 
an effective attractive interaction between the fermions, 
which is reflected by a negative sign of the scattering 
amplitude f(k) at k close to the Fermi momentum kp as 
follows from Eqs. (32 1 and (33 1. We thus employ a BCS 
approach [35] 



and reduce H to a bilinear form: 



Ubcs = J d 2 r^(r) (-^ 2 - m) ^(r) 

+ * [ d 2 rd 2 r'[A*(r,r')^(r')^(r) + h.c.}. (40) 



The order parameter (gap) obeys the gap equation 



Afr, 



F (r-r')(V>(r'Mr)}, 



(41) 



with the symbol ( ) denoting the expectation value. 
Strictly speaking, in two dimensions this approach is ac- 
curate at T = 0. At non-zero temperatures the long- 
range order is destroyed by long-wave thermal fluctua- 
tions of the phase and one has only an algebraic order. 
However, in the weakly interacting regime the character- 
istic phase coherence length is exponentially large, and 
on a distance scale smaller than this length one may still 
use the order parameter (41) and the BCS approach. In 



the uniform case this order parameter depends on r and 
r' only through the difference (r — r'). 

The Hamiltonian %bcs is reduced to a diagonal form 
Tines = Sk e *^k^k + const, by using the Bogoliubov 
transformation 



[uk exp(zkr)& k + u£ exp(— ikrW]. 



Here S is the surface area of the system, while 6k, &k 
and e k = + |Ak| 2 are annihilation/creation opera- 
tors and energies of single-particle excitations. The func- 
tions «k and t>k satisfy the well-known Bogoliubov-de 
Gennes equations. They are given by 



\/2efc(£fc + e k ) 



y/2e k fe + e k y 



where Ak = J d 3 rA(r — r')exp[ik(r — r')] is the 
momentum-space order parameter, = H 2 k 2 /2m — \i, 



and the chemical potential in the weakly interacting 
regime is close to the Fermi energy Ep — h 2 k 2 7 /2m. In 
the momentum space, the gap equation (41) takes the 
form: 



A k = - 



d 2 q 

(2tt) 2 



Vb(q-k)A q ^(g), 



(42) 



where Vo(q) is the Fourier transform of the potential 
V (r), and K{q) = tanh(e g /2r)/2e g . 



We now renormalize the gap equation ( 42 ) by express- 



ing the Fourier transform of Vo(r) through the off-shell 
scattering amplitude. The relation between these two 
quantities reads [34] : 

M 1 V ' J (2vr) 2 2(E k -E q -iO) V ' 



Multiplying Eq. (43) by JC(k')A-k> and integrating over 
d 2 k', with the help of Eq. (42) we then obtain: 



A k = -|/( k ' k ) A k' 



1 

2{E k ,-E k -iG) 



d 2 k' 



(44) 



Note that in contrast to the commonly used renor- 
malization procedure expressing the Fourier transform 
of the interaction potential through the vertex function 
r(k, k',_E) [37] at an arbitrarily chosen energy E, which 
coincides with /(k',k) for E = E k >, we use the off- 
shell scattering amplitude from the very beginning of the 
renormalization. 

At T = we put tanh(e fe /2T) = 1 and, hence, 
K,(k) = l/2e k . We then perform an analysis assuming 
that in the weak coupling limit the main contribution to 
the integral in Eq. ( 44 ) comes from momenta k! close to 
kp. It shows that the dominant pairing instability is in 
the channel with orbital angular momentum 1 = 1, since 
for higher angular momenta the interaction (scattering) 
amplitude is much smaller. The most stable low temper- 
ature phase has p x ± ip y symmetry, following from the 
fact that this phase fully gaps the Fermi surface, in con- 
trast to competing phases [35] . A full numerical solution 
of the regularized gap equation confirms this analysis. 

In fact, equation (|42|) and, hence, equation (Eil are 



not sufficient for obtaining a correct result for the order 
parameter. One should calculate the quantity 6V(q, k) 
originating from many-body effects and add it to Vb(q — 
k) in Eq. (42). The quantity (W(q, k) is a correction 
to the bare interparticle interaction Vo, and the leading 
terms of SV(q, k) are second order in Vq [3J5]. The corre- 
sponding diagrams are shown in Fig. [6] and they are the 
same as in the case of superfiuid pairing between identical 
dipolar fermions in three dimensions [40] . They describe 
processes in which one of the two colliding particles po- 
larizes the medium by creating a particle-hole pair. In 
Fig. [6]ja) the particle-hole pair then annihilates due to 
the interaction with the other colliding particle. In[6^b), 
|6][c), and[6]l the hole annihilates together with one of the 
colliding particles. In[6|b) and[6jc) the particle-hole pair 
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FIG. 6: The lowest order many-body corrections to the ef- 
fective interparticle interaction. 



is created due to the interaction of the medium with one 
of the colliding particles, and the hole annihilates with 
the other colliding partner. In |6|d) these creation and 
annihilation processes involve one and the same colliding 
particle. Including these many-body effects, which were 
introduced by Gor'kov and Melik-Barkhudarov [39], the 
gap equation becomes: 



At 



-|/(k',k)A k , 



IC(k') 



1 



SV(k',k)fC{k')A k> 



2{E k , 

d 2 k' 

(2tt)2- 



E k - iO) 



d 2 k! 
(45) 



It turns out that we should also take into account the 
difference between particles and quasiparticles (single- 
particle excitations). The latter are characterized by 
the effective mass m* , and the density of states near the 
Fermi surface is 



v F 



2nh 2 



(46) 



Thus, dealing with this quantity in the gap equation we 
should replace m with m*. Note that for short-range 
potentials the effective mass correction can be neglected. 
However, for fairly long-range interactions like l/r 3 in 
2D, this correction becomes important as was recently 
demonstrated for the case of s-wave pairing in bilayer 
dipolar systems [4"T] . 

According to the Landau Fermi-liquid theory, the 
derivative of the quasiparticle energy e(fc) with respect 
to momentum is given by |35j : 



de(k) 
dk 



h 2 k 



J F(k,k>) 



,.dN(k') d 2 k' 



5k' (2tt) 2 



(47) 



where N{k') is the distribution function of quasiparticles 
for which we take a step function N(k') = 9(kp — k'), 
so that dN(k')/dk' = -(k'/k')5(k' - k F ). The leading 
contribution to the interaction function of quasiparticles, 
F(k, k'), comes from the scattering by the l/r 3 tail of 



the interaction potential Vo(r). This contribution can be 
calculated in the first Born approximation and, including 
all (odd) partial waves, for small momenta it is given by 

F(k,k') = [V (0)-V (k-k')} = - — |k-k'|rv (48) 

m 

This simply follows from the fact that in our weakly in- 
teracting system of identical fermionic particles the in- 
teraction energy par unit area is 

Eint = \ ^o^ 3 ~ k 2)(4 1 +k 2 -k 3 «k 3 ak 2 ak 1 ) 



\ E wo) 

ki,k a 



V (k 1 -k 2 )}N{k 1 )N(k 2 ), 



with and a being the (quasi)particle operators. 
For k near kp equation (47) then yields: 



Vp 



Tikp 
m 



2n 



F(2fc F |sin ^/2|)cos^ : 



4tt 2 



(49) 



where <j> is the angle between the vectors k and k', and 
Vp = de/hdk\k=k F is the Fermi velocity. The effective 
mass is defined as = Tikpjvp. Then, using Eq. (48) 
we have: 



(50) 



m Sir 

The Gor'kov-Melik-Barkhudarov corrections and the 
replacement of the bare mass m by m* do not change our 
conclusion that the ground state has p x ± ip y symmetry. 



VI. ORDER PARAMETER AND TRANSITION 
TEMPERATURE 

In the 2D geometry that we consider, the transition 
temperature T c of a Fermi gas from the normal to su- 
perfluid regime is set by the Kosterlitz-Thouless transi- 
tion. However, in the weak coupling limit the Kosterlitz- 
Thouless temperature is very close to T c calculated in 
the BCS approach The latter follows from Eq. d45b 



as the highest temperature at which this equation has a 
non-trivial solution for the order parameter. Using the 
Px + ipy symmetry we write Ak = A(fc) exp(i0k)- Then, 
multiplying both sides of Eq. (45) by exp(— i<j> k ) and in- 
tegrating over d(j)k and d^ we obtain the same equa- 
tion (45) in which Ak and Ak' are replaced by A(fc) and 
A(fc'), the off-shell scattering amplitude /(k',k) is re- 
placed by its p-wave part f(k', k) defined in Eq. ( |13[ ), and 
6V(k',k) by its p-wave part SV(k',k). Calculating the 
contribution of the pole in the second term in the square 
brackets and expressing the off-shell scattering amplitude 
through f(k',k) by using Eq. (17), we obtain: 



A(fc) = -P f(k',k)A(k') 



lC(k')- 



1 



SV(k',k)JC(k')A(k') 



2(Ek' —E k ) 
d 2 k' 



(2n) 2 



d 2 k! 
(51) 
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FIG. 7: (color online). The zero-temperature gap Ao(fc) in 
units of Ao(kp) as a function of k/kp for Q, = 0.25, ft = 0, 
and kpr t — 0.4 near the minimum of inelastic losses (rg/Xs = 
10.5). The solid curve is the result of full scale numerics on 



(531 



the basis of Eq. (511, and the dashed curve is the result of 
Eq. " 



|A(fc)|,T c <C Ep, we divide the region of integration over 
k' into two parts: \E^ — Ep\ < u, and \Ef-i — E F \ > u, 
where |A(fc)|,T c <C ui -C E F . In the second region 
we may put ey — \£,k'\ and tanh(efc</2T) = 1, so that 
IC(k') = 1/2|&/|. Then, making use of Eq. ((52) for A(fc') 
we see that the integration in this region gives A(k F )A, 
where the quantity A is temperature independent. Divid- 
ing both sides of the gap equation by A(k F ) we obtain: 



- 2 [y-(k F ) + sv(k F , kF)] t^( v^gggg 



27T?l Jo 



2T 



■4 = i. 



(54) 



As T — > T c the order parameter tends to zero and we 
may put A(k F ) — in Eq. (54). Then, subtracting this 
equation at T — > T c from the equation at T — we have: 



\tanh{£ k ,/2T c ) 



1 



v/e + |A (MI 2 



= o. 



The integral converges at £fc< of the order of T c or 
|Ao(fci?)|, and we may extend the upper limit of inte- 
gration to infinity. This gives: 



The symbol P denotes the pr inci pal value. The ampli- 
tude f(k',k) is real, and Eq. (511 is convenient for ana- 
lytical and numerical calculations of the order parameter 
and T c . 

The leading contribution to the integral on the right 
hand side of Eq. ( 51 ) is related to the first term in square 



brackets and comes from a narrow vicinity of the Fermi 
surface. Omitting other contributions we then establish 
a relation between A(fc) and the order parameter on the 
Fermi surface: 



A(fc) = A(fc F 



f(k F ,k) 
Rk F ) • 



(52) 



For small momenta it is sufficient to use Eq. (34) for 
f(k F ,k) and Eq. ((32) for f(k F ). This yields: 



A(fc) = 



ZnAjkp) f^F (-1,5,2, g-), k 



< k f 



f {-hh^w h k>kF 



(53) 



so that A(fc) oc k for k < k F , and it becomes k- 
independent for k > k F . In Fig. 7 we compare the result 
of Eq. ( 53 ) with the result of full scale numerics at T = 
for the interaction potential Vo(r) at = 0.25(5 and 
= 0. 



Equation ( 52 ) can be used to calculate the contribu- 



tion to the integral in Eq. (51) from momenta k! away 
from the Fermi surface. This immediately leads to a 
relation between the zero-temperature order parameter 
on the Fermi surface, Ao(fci?), and the critical tem- 
perature T c . Putting k — k F in Eq. (51) and taking 



into account that in the weak coupling regime one has 



|Ao(fcfOI, 



(55) 



which is the same relation as in the case of the s-wave 
pairing in a two-component weakly interacting 3D Fermi 
gas [3S]. Note that equation (54) allows one to establish 
a relation between T c or |Ao(Ff)| and |A(fcp)| at any 
temperature. Having in mind that the momentum de- 
pendence of the order parameter follows from Eq. ( 52 ) 



we see that the calculation of T c provides us with the 
order parameter at any k and temperature. 

We now calculate the critical temperature T c on the 
basis of the gap equation (51 ) at k = hp. In the first line 



of this equation we again divide the region of integration 
into two parts: \Ek> — Ep\ < w, and \Ek> ~ Ep\ > uj. In 
the first region we put A(fc') = A(fcp) and f(k',k F ) = 
f(k F ). The contribution of the second term in square 
brackets is then equal to zero. In the first term we set 
e k , = |£ fc ' | and, hence, JC{k') = (2|&/ 1)" 1 tanhflfr \/2T c ). 
Denoting the result of the integration of this term as 
Ai(k F ) we have: 



A 1 (kp)=A(k F ) 



4fc F r* 
3tt 



1 



3tt 
16 



kpr* ln(pk F r*) 



x In 



2e C LO 
ttT c 



(56) 



where we used f(k F ) 
given by equations (|32|) 



Mkp) + h{kp), with /i, f 2 
(|33|) and the numerical factor 



p introduced after Eq. (1311) 

In the second region, \Ey — E F \ > lu, we put JC(k') 
l/2|£fc/|. The contribution of this region is small as 
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~ fc F r* compared to the result of Eq. (56). Therefore, 
it is sufficient to retain only the leading low-momentum 
contribution to the off-shell amplitude f(k',kp). Thus, 
using Eq. Q for f(k',k F ) and Eq. Q for A(k') we 
write this contribution as 



A 2 (fc F ) = A(fc F ) 



37rr* 
8k F 



" k' 3 dk' 
Wf - k' 2 ) 



Il k 12 

2' 2' ' kj 



where k^ = y2rn(Ep — lu)/Tl 2 . Taking into account that 
lu <C Ep we then obtain: 



A 2 (fc F ) = A(k F ) 



where 
i] = 1 



97T 2 r 1 
64 Jo 

xdx 



F 2 
0.78. 



3tt 



--, -,2,x 
2' 2' ' 





- V 







(57) 



F 2 



2' 2' ' 



1 - a; 

The main contribution to the term in the second line 



of Eq. (51), describing Gor'kov-Melik-Barkhudarov cor- 
rections, comes from the vicinity of the Fermi surface, 
i.e. from the region where \Ey — Ep\ < u). The result of 
the integration reads: 



-^(M o- ln -^f- ■ ( 58 ) 



2tt 



In the Appendix we show that in the low-momentum 
limit one has 



SV(k F ,k F ) = a — (fc F r*) 
m 



(59) 



where a ~ 2.3. 



Then, making a summation of the contributions (56), 



we obtain 

4fc F r* 
3vr 



(57), and (58), and dividing the gap equation by A[kp 
2e c £ F 



In 



(fc F r*) 2 ln 



Ep 



{ l -Hpk F r,) + ^- 



= 1. 



(60) 



Note that we put uj ~ -E F in the argument of the loga- 
rithm in the Gor'kov-Melik-Barkhudarov term and in the 
term proportional to (fc F r + ) 2 ln(pfc F r*). This is justified 
because ui can be chosen as a small numerical fraction of 
Ep, and these terms are small as ~ kpr* compared to 
the leading term. 

We now recall that the bare mass m should be re- 



placed by the effective mass following from Eq. ( 50 1 . 
The relative difference between and m is propor- 
tional to kpr* and is small. Therefore it is sufficient 
to replace m with m* only in the multiple r, a m in 



the first term of Eq. (60). Using Eq. (50) this leads to 



the appearance of a new term (16/97r 2 )(fc F r st ) 2 ln(Ep/T c ) 
in equation (60), which is equivalent to replacing a by 
a = a — 32/97T ~ 1.17. Relying on the inequality 
fcpr, < 1 we then immediately find the transition tem- 
perature 



E F - (fc F r,)9- 2 /64 ex P ( 4 fc F rJ ' 



where 



9ir 2 9ir (2 

'' < X1>< ' ~~M^ np ~ 32 5 + C ~ V + ln 



Substituting the values of p, 77, and a specified above we 
obtain: 



K ~ 0.16 exp 



9n 3 A 
64 



(62) 



Let us make several important statements regarding 
equation (61). First of all, the factor in the exponent is 
~ — (/c F r*) _1 , and the ratio T c /Ep can be made larger 
than 10 -2 even for relatively small fc F rv This is a con- 
sequence of the anomalous scattering due to the attrac- 
tive 1/r 3 tail of the potential Vq (r) , which gives a scat- 
tering amplitude proportional to k. Moreover, an ac- 
curate calculation of the scattering amplitude reveals a 
(kr*) 2 In(Ax*) term, which leads to the appearance of a 
power law factor (£; F r >t )~ 97r Z 64 in front of the exponent 
and provides an additional increase of the critical temper- 
ature. This behavior is in contrast to the p-wave pairing 
of short-range interacting atoms, where the factor in the 
exponent is inversely proportional to —kp and the BCS 
critical temperature is vanishingly low. 

The dependence of T c on the short-range behavior of 
the potential Vo(r) is contained in the factor k through 
the coefficient A. In Fig. 8 we present A and k as func- 
tions of rs/Xs for £Ir — 0.25(5 and f3 ranging from to 
0.2. These results show that K can be varied by two or- 
ders of magnitude by changing rs/Xs [i.e. by changing 
the depth of the potential Vo(r)] while remaining in the 
regime without bound states in this potential, such that 
a rapid three-body decay is absent. 

In Fig. 9 we display the ratio T c /Ep as a function of 
kpr* for f2 F = 0.255 and rs/Xs close to the minimum 
of the loss rate at a given f3. We compare the full scale 
numerical solution of the gap equation (51 ) with the ana- 



lytic expression, Eq. ( 60 1 . The discrepancy between ana- 
lytic and numerical results as &; F r* approaches 1 is to be 
expected as the analytic method is perturbative in this 
parameter. From the results of Fig. 9 we conclude that 
a critical temperature of the order of 5% of the Fermi 
energy is realistic. For typical 2D densities n ~ 10 8 — 10 9 
enr 2 the Fermi energy of alkali atom molecules is of the 
order of hundreds of nanokclvins. Then, as we see from 
Fig. 9, the superfluid transition temperature can be as 
high as 10 or 20 nK. 
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FIG. 8: (color online). The coefficients A and k versus rs/Xs 
for fl R /S = 0.25 and (3 = (black, solid), 0.1 (blue, dashed), 
and 0.2 (red, dash-dotted). The coefficients are presented for 
the values of rs/Xs at which bound states in the potential 
Vo(r) do not exist. 
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FIG. 9: (color online). Critical temperature T c in units of the 
Fermi energy as a function of kpr» for Qr/8 = 0.25 near the 
minimum of inelastic losses. In (a) /3 = 0, rs/Xs ~ 10.7, in 
(b) P = 0.1, r 6 /\$ = 11.4, and in (c) /3 = 0.2, r s /X s = 13.2. 
The data points are obtained from full scale numerics on the 
basis of Eq. ( |51[ ), and the curves from Eq. ( |60[ ) 43 . 



VII. CONCLUDING REMARKS 

The results described in this paper indicate that it is 
realistic to create the superfluid topological p x +ip y phase 
with alkali atom polar molecules. The key point is the 
anomalous scattering by the attractive 1 /r 3 tail of the in- 
teraction potential Vo(r), which leads to a relatively large 
(and negative) p-wave scattering amplitude and thus to 
an achievable transition temperature T c even in the BCS 
regime. What we have shown in the present paper is the 
possibility to manipulate T c by modifying the short-range 
part of the potential Vo(r). As is seen from Fig. 8, the 
pre-exponential coefficient k can be varied within 2 or 3 
orders of magnitude, and we may still remain in the BCS 
regime with fairly small inelastic losses. 

Another result of this paper is that the addition of a 
dc electric field reduces the inelastic losses, and one can 
work at higher densities. Consider, for example, 6 LiCs 
molecules which have a permanent dipole moment of 5.5 
D [441 05]. For n R = 0.12(5 and /3 = 0.165 the dipole- 
dipole distance is r% = 100 nm and kpr* = 1 at 2D 
densities of 10 9 cm -2 corresponding to a Fermi energy 
of 200 nK. Then, our calculations give the rate constant 
10~ 9 cm 2 /s for rg/Xg — 39, which for the selected £Ir/5 
and (3 corresponds to the absence of bound states in the 
potential Vo(r) and is not far from the minimum of in- 
elastic losses. Thus, the lifetime of the gas is about 1 
s at these densities. With kpr* = 1 we are beyond the 
BCS regime. However, the results of Fig. 9 indicate that 
in this case we are likely to have a critical temperature 
of at least a twentieth part of Ep, which is T c ~ 10 nK. 
Note that the ratio rg/Xg = 39 corresponds to rg ~ 15 
nm and using the rotational constant B ~ 0.27 K [?5] we 
find that j3 = 0.165 corresponds to Ed c — 350 V/cm. 

In order to provide the 2D regime of scattering and su- 
perfluid pairing one should have the confinement length 
in the ^-direction, l z < rg. For the above example of 
6 LiCs molecules with rg ~ 15 nm this requires a con- 
finement frequency ui z — fi/ml 2 z of the order of 300 kHz. 
Such very high frequencies are not unrealistic and are 
achievable by a tight optical confinement. In particular, 
frequencies close to 100 kHz have been used for sideband 
cooling of cesium atoms [4BTl48j . 

Similar estimates are obtained for LiK and NaK 
molecules which have dipole moments of 3.5 and 2.7 D, 
respectively [H]. Here one can have the Fermi energy 
above 300 nK (n > 10 9 cm -2 ) and T c on the level of 
tens of nanokelvins, with a lifetime of the order of a sec- 
ond. However, the necessary confinement frequencies are 
about 500 kHz. 

The requirement for the confinement frequency can be 
relaxed by introducing a shallow optical lattice and thus 
increasing the effective mass of the molecules. For exam- 
ple, increasing the effective mass by a factor of 5 allows 
one to increase rg to 25 nm for 6 LiCs molecules, which 
corresponds to the confinement frequency of 100 kHz. 
Then, however, we have to have larger r + and even for 
kpr* ~ 1 the Fermi energy is about 20 nK, so that T c 
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will be on the level of nanokelvins. Similar estimates are 
obtained for LiK and NaK molecules where under the 
same increase of the effective mass the confinement fre- 
quency can be decreased to 200 kHz, with T c of the order 
of nanokelvins. Note that increasing the effective mass 
by a factor of 10 makes the situation promising even for 
40 KRb molecules which have a much smaller dipole mo- 
ment d = 0.6 D pQ. Then, taking r$ ~ 20 nm, which cor- 
responds to the confinement frequency of 200 kHz and 
to = 100 nm we find T c on the level of nanokelvins 
(Ep ~ 30 nK) at densities of 10 9 cm~ 2 and with a life- 
time of the order of seconds. 

Note that the strong confinement in the z-dircction 
with frequencies lo z exceeding 100 kHz, originating from 
the condition l z % r$, is only necessary for consider- 
ing the relative motion of colliding molecules as two- 
dimensional at intermolecular distances r ~ rg and re- 
vealing the role of short-range physics for the elastic p- 
wave interaction and inelastic processes. The 2D super- 
fluid pairing by itself, which leads to the p x + ip y symme- 
try of the ground state and mostly originates from the 
anomalous scattering, simply requires the dipole-dipolc 
distance 7% significantly larger than l z . For all examples 
given above, this condition is well satisfied at a frequency 
lu z ss 50 kHz. However, the absence of the inequality 
fs ~ h. introduces a numerical uncertainty in the coef- 
ficient k in Eq. (61) and can make the inelastic decay 
faster. 

The formation of the p x + ip y superfluid phase should 
be apparent in numerous observables. Superfluidity it- 
self can be detected by the means that have been used to 
detect s-wave superfluids [2PJ ISO] • The most interesting 
new properties arise in the presence of quantized vor- 
tices, which are predicted to carry zero energy Majorana 
modes on their cores [19] . Vortices can be induced in su- 
perfluid Fermi gases by causing the gas to rotate [I5J OjT] • 
The appearance of zero energy Majorana modes on the 
vortex cores will lead to signatures in the RF absorption 
spectrum: the Majorana modes are predicted to lead to 
a series of sharp lines below the superfluid gap. These 
lines arise from processes in which a fermion is excited 
out of the Majorana modes [52i. Ultimately one would 
hope to probe non-abelian exchange statistics of these 
vortices 1531. 
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Appendix 

In this Appendix, we demonstrate how the Gor'kov- 
Melik-Barkhudarov corrections to the bare interparticle 
interaction may be computed. The corrections consist of 
the four diagrams depicted in Fig. 6 [ID]. In the low- 
momentum limit they are dominated by the long-range 
1/r 3 tail of the potential Vo(r) and, hence, we may ap- 
proximate this potential as 

Mr) « { °l { tf /rn)rm/l * t T r t? ' (63) 

where ro is much smaller than the de Broglie wavelength 
of particles. The Fourier transform of the potential takes 
the form: 



Vo(q) 



2irh 2 r* 2nh 2 , 



qr. 



(64) 



to ro m 
The expressions for the four diagrams of Fig. 6 are 

d 2 q 7V(q + p_/2)-7V(q-p_/2) 



SV a (p',p) 



SV b (p',p) 



(2tt)= 





£q+p_/2 - £q-p_/2 



xV 2 (p 



(65) 



d 2 q jV(q + p_/2)-jV(q-p_/2) 

( 27r ) 2 Cq+p_/2 - Cq-p_/2 

xU (p_)F (q-p+/2), (66) 



d 2 q A(q + p_/2)-JV(q-p_/2) 



( 27r ) 2 £q+P-/2 - £q-p_/2 

xV (p_)Vo(q + P+/2), (67) 



W,P) = -/^ 



d 2 q 7V(q + P+ /2)-A(q-p + /2) 



£q+p + /2 - £q-p+/2 



xV„(q-P-/2)Vo(q + P-/2), 



(68) 



with p± = p' ± p, and £ q = h q 2 /2m — fi. In the weakly 
interacting regime we may put the chemical potential fi 
equal to the Fermi energy, and at very low temperatures 
T < T c we may take the step function N(q) = Q(kp — \q\) 
for the Fcrmi-Dirac distribution. 



The p-wave part of the diagrams (|65|-(68| is 



SV(p',p)=J28V J (p', P )=r^e-^ ]r^(p',p), (69) 

where <f> is the angle between p' and p. One can eas- 
ily check that the momentum- independent term of Vq (q) 
does not contribute to the sum of the four diagrams. 

For the analytical calculation of T c we only need the 
contributions SVj(p' 1 p) on the Fermi surface, i.e. for 
p' = p = kp. In this case each contribution can be repre- 
sented in the form SVj(kp, kp) = ctj(h 2 /m)(kpr*) 2 . The 
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calculation of the diagram |6ja) is straightforward and it 
gives a a = 2tt. The other diagrams we calculated numer- 
ically. The values of the coefficients are 



a a = 2tt, at — a c ~ -1.5, 



ad 



-1.0, 



and thus we find the result of Eq. (59): 



h 

SV(k F ,k F ) = a — (k F r*) , 
m 



(70) 



(71) 



with a ~ 2.3. 

For the numerical solution of the gap equation we 
used a simplified dependence of 5V on p' and p, namely 



the dependence following from the diagram (a) with 

Vb (| P | ) = (27r7i /m)|p_|r*. This dependence is given 

by: 



SV{p',p)=5V(k F ,k F ){^ + 



d4 
27 



0(p_ — 2kp) cose 



P-Jp 2 _~4k 2 F 



}. (72) 



For (p+p') < 2k F we have SV(p',p) = SV(k F ,k F )pp' /k F , 
and this quantity rapidly decays at larger momenta. 
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